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ABSTRACT 



Aims. Nova Cyg 2006 has been intensively observed throughout its full outburst. We investigate the energetics and evolution of 
the central source and of the expanding ejecta, their chemical abundances and ionization structure, and the formation of dust. 
Methods. We recorded low, medium, and/or high-resolution spectra (calibrated into accurate absolute fluxes) on 39 nights, along 
with 2353 photometric UBVRcIc measures on 313 nights, and complemented them with IR data from the literature. 
Results. The nova displayed initially the normal photometric and spectroscopic evolution of a fast nova of the Fell-type. Pre- 
maximum, principal, diffuse-enhanced, and Orion absorption systems developed in a normal way. After the initial outburst, the 
nova progressively slowed its fading pace until the decline reversed and a second maximum was reached (eight months later), 
accompanied by large spectroscopic changes. Following the rapid decline from second maximum, the nova finally entered the 
nebular phase and formed optically thin dust. We performed a photo- ionization analysis of the emission- line spectrum during 
the nebular phase, which showed a strong enrichment of the ejecta in nitrogen and oxygen, and none in neon, in agreement with 
theoretical predictions for the estimated 1.0 Mq white dwarf in Nova Cyg 2006. The similarities with the poorly investigated 
V1493 Nova Aql 1999a are discussed. 

Key words. Stars: novae - Stars: individual: V2362 Cyg 



1. Introduction 

Classical novae are powered by thermonuclear runaways 
occurring on white dwarfs in close binary systems of the 
cataclysmic variable type, during which H-rich material 
accreted from a low-mass companion is burned and ex- 
plosively ejected into the interstellar medium, typically at 
^1000 km sec -1 bulk velocity and 10 -4 M Q total mass. 
Mixing with the underlying white dwarf material affects 
the chemistry of the ejecta, which is dominated by non- 
equilibrium CNO-burning products. Classical novae are 
instrinsically very luminous objects (with a luminosity ap- 
proaching or exceeding the Eddington luminosity for a 
1 M Q object, and My staying brighter than —5.5 mag 
for the first two weeks past maximum) , and are easily ob- 
served throughout a vast fraction of our Galaxy, the Local 
Group, and beyond it up to the Virgo Cluster. The well- 



* Tables 1, 3 and Figures 7, 8, 16-20 available only electron- 
ically 

** Based in part on observations obtained with the Asiago 
telescopes 



established relation between their absolute magnitude and 
the speed of decline makes novae useful distance indicators 
(e.g. Cohen 1985). 



Significant research carried out in the past at the 
time of photographic spectroscopy to qualitatively char- 
acterize the spectroscopic evolution of novae (e.g. Payne- 
Gaposchkin 1957, hereafter PG57, and references therein) 
has led to the discovery of common paths followed by the 
majority of novae, irrespective of their speed class and par- 
entship with galactic stellar populations (eg. McLaughlin 
1960, hereafter McL60). Analogous and more quantita- 
tive investigations in the era of CCD detectors are ap- 
preciably rarer due to the large amount of telescope 
time allocation required, expecially if observations aim to 
cover the entire nova outburst. Similarly, historical photo- 
graphic or amateur visual lightcurves exist for many novae 
(e.g. Mobberley 1999, Kiyota 2004), but densely mapped, 
multi-band, and accurate CCD lightcurves extending from 
maximum to the faintest evolutionary stages are quite 
rare. 
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Table 1. Portion of Table 1, available only electronically, 
to illustrate its content. The quoted error is the Poisson 
error and the last column provides identification of the 
telescope (see details in Sect. 2.1). 



HJD 



band 



mag 



id 



2454020 


2807 


B 


11 


784 


0.002 


d 


2454020 


2850 


U 


11 


891 


0.005 


e 


2454020 


2887 


V 


11 


611 


0.002 


d 


2454020 


2933 


Rc 


10 


785 


0.001 


d 


2454020 


2977 


Ic 


10 


660 


0.007 


e 


2454020 


2978 


Ic 


10 


634 


0.001 


d 


2454020 


3067 


V 


11 


644 


0.006 


e 


2454020 


3178 


B 


11 


869 


0.007 


e 


2454020 


3253 


Rc 


10 


912 


0.004 


e 


2454020 


3470 


U 


11 


741 


0.028 


b 


2454020 


3539 


B 


11 


909 


0.022 


b 


2454020 


3574 


V 


11 


633 


0.022 


b 


2454020 


3577 


Rc 


10 


817 


0.012 


b 


2454020 


3580 


Ic 


10 


606 


0.032 


b 


2454021 


2753 


Rc 


10 


832 


0.004 


e 


2454021 


2820 


B 


11 


751 


0.006 


e 


2454021 


2918 


V 


11 


570 


0.005 


e 


2454021 


3003 


Ic 


10 


591 


0.007 


e 


2454021 


3102 


B 


11 


723 


0.002 


d 


2454021 


3179 


V 


11 


562 


0.002 


d 


2454021 


3183 


u 


11 


760 


0.038 


e 


2454021 


3231 


Rc 


10 


745 


0.001 


d 


2454021 


3283 


Ic 


10 


593 


0.001 


d 


2454021 


3755 


U 


11 


605 


0.026 


b 


2454021 


3824 


B 


11 


812 


0.014 


b 


2454021 


3859 


V 


11 


566 


0.030 


b 


2454021 


3862 


Rc 


10 


795 


0.020 


b 


2454021 


3865 


Ic 


10 


556 


0.015 


b 



This paper presents new and extensive multiwave- 
length photometric and spectroscopic observations of 
Nova Cyg 2006 (= V2362 Cyg, hereafter NCyg06). It was 
discovered at magnitude 10.5 by H. Nishimura on photo- 
graphic plates exposed on April 2.807 UT (as reported by 
Nakano 2006), while the nova was still on the rise to its 
maximum. Within a few hours of the discovery announce- 
ment, several spectroscopic confirmations of the nature of 
the nova were made, as summarized by Yamaoka (2006a). 
This nova is characterized by one of the most unusual pho- 
tometric behaviors ever recorded for a classical nova. We 
investigate in detail the initial outburst, the rebrighten- 
ing that occurred several months later, and the advanced 
decline phase. We also report photo-ionization modeling 
of the emission spectrum from the expanding ejecta to es- 
timate their total mass, their chemical partition, and the 
temperature and luminosity of the ionizing central white 
dwarf, thus adding to the limited number of novae (~30) 
for which such analysis has been performed. 
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Fig. 1. Lower six panels: photometric evolution of Nova 
Cyg 2006 in the U BV RqIc bands. Upper six panels: tem- 
poral evolution of [OI] 8446/Ha and [OI] 6300/ [OI] 6364 
flux ratios, of [OI] 6300, [Oil] 7325, [OIII] 5007, and Ha 
integrated line flux (expressed in magnitudes relative their 
peak values 9.364xl0~ 12 , 7.151xl0~ 12 , 1.249X10" 11 , and 
3.351xl0~ 10 erg cm~ 2 sec -1 , respectively), of Ha equiv- 
alent width, and of FWHM (in km/s) of H/3, [OI] 6300, 
and [Nil] 5750 A emission line profiles. After day +250, a 
growing fraction of [Nil] 6548, 6584 emission contributed 
to the Ha flux and equivalent width. 
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2. Observations 

2.1. Photometry 

We collected (cf. Table 1) CCD and photoelectric pho- 
tometry of NCyg06 with a number of telescopes, densely 
covering all stages of its temporal evolution. All obser- 
vations were reduced against the UBVRcIc comparison 
sequence presented by Frigo et al. (2006). This sequence 
spans a wide range in color and allows suitable transforma- 
tion from the local to the standard systems. The different 
instruments are associated in the next paragraph to let- 
ters a to g, which arc used in Table 1 and throughout the 
rest of the paper to identify them. 

Observations were obtained with: (a) the Sonoita 
Research Observatory (SRO) 0.35-m Celestron C14 
robotic telescope using BVRqIc Optec filters and an 
SBIG STL-1001E CCD camera, 1024x1024 array, 24 
/iin pixels = 1.25"/pix, with a field of view of 20'x20'; 

(b) the 0.42-m f/5.4 Newtonian telescope operated 
by Associazione Ravennate Astrofili Rhcyta in Bastia 
(Ravenna, Italy), equipped with an Apogee Alta 260e 
CCD camera, 512x512 array, 20 /im pixels = 1.83"/pix, 
field of view of 16'xl6' and Schuler UBVRcIc filters; 

(c) the 0.50-m f/8 Ritchcy-Chretien telescope operated on 
top of Mt. Zugna by Museo Civico di Rovereto (Trento, 
Italy) and equipped with Optec UBVRcIc filters. The 
CCD is an Apogee Alta U42 2048x2048 array, 13.5 
pixels = 0.70"/pix, with a field of view of 24'x24'; (d) 
the 0.30-m Meade RCX-400 f/8 Schmidt-Cassegrain tele- 
scope owned by the Associazione Astrofili Vallc di Cembra 
(Trento, Italy). The CCD is an SBIG ST-9, 512x512 array, 
20 /jm pixels = 1.72"/pix, with a field of view of 13'xl3'. 
The B filter is from Omega and the V^i?c-^c filters from 
Custom Scientific; (e) the 0.13-m f/6.6 Vixen ED130SS re- 
fractor privately owned by one of us (G.C.) and operated 
in Trieste (Italy). It is equipped with Custom Scientific 
UBVRcIc filters and a Starlight SXV-H9 CCD camera, 
1392x1040 array, 6.45 ^m pixel = 1.55"/pix for a field 
of view of 36'x27'; (/) the 0.28-m Celestron Cll tele- 
scope privately owned by one of us (S.D.) and operated 
in Cembra (Trento, Italy). It is used with an Optec SSP- 
5 photoelectric photometer and standard B,V Johnson 
filters; and finally (g) the 1.0-m Ritchey-Chretien tele- 
scope of the U. S. Naval Observatory, Flagstaff Station 
(NOFS). A Tektronix/SITe 2048x2048 thinned, backside- 
illuminated CCD was used. The telescope scale is 0.6763 
arcsec/pixel, with a field of view of 23'. 1x23'. 1 arcmin. 

The 2353 photometric measures (89 in U, 614 in B, 621 
in V, 498 in R c , and 531 in I c ) we collected on NCyg06 
on 313 different nights arc listed in Table 1 (available elec- 
tronic only). The overall nova light-curve and color evolu- 
tion are shown in Fig. 1 (see also Kimeswenger et al. 2008), 
using data only from instruments a, d, and g (Figs. 16-20, 
available electronic only, present the equivalent plots using 
all of the data in Table 1). 



Table 2. Journal of spectroscopic observations. 



date 


UT 


HJD 


range 


A/pix 


expt 


g 


At 


AV 


ZUUU Uo 


01'44 


'iR'i'i ^70 
oooo.u / u 


mnn 7ds k ; 


1 76 


900 s 




2 5 


73 




02'56 


OOOO.UZU 


fii 7^ f,Rin 

Ull u-uoou 


32 


1200 s 


2 


2 5 


73 


04 12 


02'07 


■iR37 ^Rfi 


IQnn 797D 


1 76 


1800 s 




6 5 


1 64 




03:09 


3837.629 


6200-6800 


0.32 




2 


6.5 


1.64 


04 13 


00:21 


3838.512 


3600-7300 


22000 


1600 s 




7.4 


1.79 


04 14 


00:25 


3839.515 


3600-7300 


22000 


1600 s 


$ 


8.4 


1.91 


04 24 


00:54 


3R4Q 


l^n^ 7780 
OOUU / / ou 


4.24 


70 s 


Q 


18.4 


2.70 




01:00 


'iRAQ ^4n 


filQ^i 7n4 K l 


0.64 


150 s 


g 


18.4 


2.70 




01:06 


3849.544 


S97 1 ! Q 1 70 
o^> i u-y i / u 


0.89 


150 s 


Q 


18.4 


2.70 


04 27 


01:54 


3R^9 ^77 


i^n^ 77sn 

OOUU / / ou 


4.23 


120 s 


g 


21.5 


2.89 




02:04 


^Rt9 ^RA 


filQ 1 ! 704^ 


0.64 


240 s 


g 


21.5 


2.89 




02: 18 


^Rt9 ^Q4 


S97t Q 1 70 
o^> i u-y i / u 


0.89 


420 s 


Q 


21.5 


2.89 


05 1 1 


23' 12 


oou / .1UU 


iR^n 79nn 


1 76 


5400 s 




36 4 


3 44 


05 25 


22' 31 


OOO 1 .^±0 


iRnn 7nnn 

oouu- i uuu 


1 76 


6300'" 




50 3 


3 74 


06 01 


02:01 


3887.584 


3495-7805 


4.24 


210 s 




56.5 


3.83 




02'11 


^RS7 
ooo / .oy i 


filQfl 704^ 


64 


180 s 


g 


56 5 


3 83 


06 15 


21'04 


^Qfl9 178 


iRnn 7nnn 

OOUU- 1 uuu 


1 76 


1800 s 




71 3 


4 00 


06 22 


21:59 


3909.417 


3800-7075 


3.52 


5400 s 


3 


78.3 


4.06 


06 24 


02:14 


3910.595 


3505-7810 


4.24 


160 s 


Q 


79.5 


4.06 




02:23 


3910.600 


6395-7045 


0.64 


270 s 


g 


79.5 


4.06 


07 10 


21:28 


QQ97 QGfi 


IR^n 71 in 

OOOU / 1 ou 


3.52 


5400 s 


3 


96.3 


4 . 1 4 


07 19 


00:20 


^ ^17 


ir,nn 7inn 

OUUU 1 ouu 


9 9 fi (1(1 


6300 s 


g 


104.2 


4.14 


07 29 


21:07 


1Q4P, 1R9 


iR^n 7i sn 

OOOU- 1 loll 


3.52 


5400 s 


3 


1 1 5.3 


4.13 




23 '09 


'iQAp, AP.7 


ijyou-uu / u 


64 


5400 s 


4 


1 1 5 i 


4 13 


08 21 


92:28 


IQfiQ 41Q 
oyuy.^oy 


iR^n r.QQn 

ooou uy y u 


3.52 


5400 s 


3 


138.3 


4.04 


08 31 


21:21 


^Q7Q 'iQ'l 
oy / y.oyo 


finnn fififi^ 

uuuu-uuuo 


0.64 


5400'' 


i 
4 


148.3 


3.99 


09 05 


21'26 


'iQRA IQfi 


irio 7nnn 

ooou- / UUU 


3 52 


5400 s 


3 


153 3 


3 96 


09 20 


20'52 


oyyy.o i o 


IR^n 7D4D 

OOOU- / U4U 


3 52 


5400 s 


3 


168 3 


3 83 


09 26 


20:53 


4005.373 


3505-7805 


4.24 


300 s 


Q 


174.3 


3.77 




21:11 


4005.386 


6395-7045 


0.64 


540 s 


g 


174.3 


3.77 


10 08 


18'50 


4m 7 9R7 


iR^n 7f\R^ 
ooou- /UOU 


3 52 


3600 s 


3 


186 2 


3 62 


10 27 


20:54 


dlTlfi 171 


iR^n 7i fin 

ooou / 1 uu 


1.76 


3000'' 




205.3 


3.34 


11 05 


20:28 


4045.354 


6130-6790 


0.32 


3600 s 


2 


214.3 


3.09 


11 12 


18:42 


4052.280 


3850-7085 


1.76 


3600 s 


I 


221.2 


2.85 


11 13 


20:19 


4n^1 9Sfl 


innn RRnn 

OUUU oouu 


3000 


3600 s 


j 


222.2 


2.81 


11 22 


20:02 


4nfi9 11^ 


1R4n 71 fin 

OdU / 1 UU 


1.76 


3600 s 


j 


231.2 


2.32 


12 10 


17:33 


4080.231 


6190-6840 


0.64 


3600 s 


i 

4 


249.1 


3.78 




18:33 


4080.273 


3850-7055 


3.52 


1800 s 


3 


249.2 


3.82 


12 13 


20'25 


4nsi i^n 


iR^n 7nT K ; 

OOOU- 1 uoo 


3 52 


2700 s 


3 


252 3 


4 79 


12 14 


17:51 


4084.243 


3505-7580 


4.24 


630 s 


6 


253.1 


4.94 




19:01 


4084.293 


3600-8900 


3000 


3600 s 


7 


253.2 


4.95 


12 21 


18:40 


4091.277 


3900-7300 


3.52 


3600 s 


3 


260.2 


5.33 




19:59 


4091.332 


6140-6800 


0.64 


2700 s 


4 


260.2 


5.33 


12 31 


17:43 


4101.237 


4000-7240 


3.52 


3600 s 


3 


270.1 


5.60 


2007 01 11 


18:50 


4112.283 


3870-7240 


3.52 


7200 s 


3 


281.2 


5.90 


01 12 


18:30 


4113.269 


3870-7240 


3.52 


4500 s 


3 


282.2 


5.91 


01 20 


17:24 


4121.223 


3505-7800 


4.24 


690 s 


6 


290.1 


5.97 




18:27 


4121.266 


3600-8900 


3000 


1800 s 


7 


290.2 


5.97 


01 28 


18:35 


4129.272 


4550-5200 


0.83 


5400 s 


10 


298.2 


6.04 


01 30 


19:17 


4131.301 


4550-5200 


0.83 


7200 s 


10 


300.2 


6.05 


03 11 


04:36 


4170.683 


3505-7800 


4.24 


1000 s 


6 


339.6 


6.35 


07 20 


00:30 


4301.523 


3505-7800 


4.24 


6000 s 


6 


470.4 


7.19 




01:48 


4301.577 


6395-7045 


0.64 


1800 s 


8 


470.5 


7.20 


09 14 


21:50 


4358.410 


3505-7800 


4.24 


7200 s 


6 


527.3 


7.01 




23:30 


4358.479 


6395-7045 


0.64 


3600 s 


8 


527.4 


7.01 


09 15 


00:17 


4358.512 


8275-9170 


0.89 


1800 s 


9 


527.4 


7.01 



HJD: heliocentric Julian day (-2450000) 
S: instrument identification code (detailed in Sect. 2.2) 
At: time elapsed since maximum V brightness (April 5.6, 2006) 
AV: decline in magnitude from maximum (V max — 8.0) 
Fifth column: spectral dispersion, except where a slanted number is given, 
which indicates a resolving power (for Echelle instrumental configurations) 



2.2. Spectroscopy 

A journal of the spectroscopic observations we ob- 
tained on NCyg06 is given in Table 2, where a num- 
ber from 1 to 10 in column seven identifies the differ- 
ent telescope/spectrograph configurations, described be- 
low. Reduction of all spectroscopic observations was car- 
ried out in a homogeneous way with IRAF, and involved 
corrections for bias, dark, and flat fields. All spectra (with 
the exception of the Echelle ones) were calibrated into 
absolute fluxes via observations on the same nights and 
instrumental configurations of one or more of the follow- 
ing spectrophotometric standard stars: HR 718, HR 5501, 
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7001, 8634, and HD 161817. On several nights, the non- 
photometric conditions required a shift be applied to the 
zeropoint of the fluxed spectra of the nova. This shift was 
determined by comparing the fluxes integrated over the 
photometric bands on the nova spectra with the corre- 
sponding fluxes from CCD photometry of the nova for the 
same night (or extrapolated from adjacent nights in one 
case). Higher resolution observations of the nova (codes 
2, 4, 5, 7, 8, 9, and 10 in Table 2) were always accom- 
panied by identical observations on the same night of one 
or more of the following stars: HD 196740, HD 209833, 
HD 212571. These are fast-rotating B stars with no Ha 
emission, and their flat, featureless continua provided an 
excellent reference background to record interfering ab- 
sorptions by telluric H 2 and 02- 

Low and medium resolution, absolutely fluxed spec- 
tra of NCyg06 were obtained with the 0.6-m telescope of 
Osservatorio Astronomico G. Schiaparelli (Varese, Italy), 
equipped with a grating spectrograph and an SBIG ST- 
10XME CCD, 2184x1472 array, 6.8 jum pixel. The slit was 
always aligned east-west, and its projected width was kept 
constant to 2.0" through the whole observing campaign. 
The spectrograph spatial scale is 1.0 "/pix. Five differ- 
ent set-ups were used: 600 ln/mm grating, range 3900- 
7100 A, lxl binning, and 1.76 A/pix (code 1 in Table 2); 
600 ln/mm grating, range 3900-7100 A, 2x2 binning, and 
3.52 A/pix (code 3); 1800 ln/mm grating, range 6200- 
6900 A, 1 x 1 binning and 0.32 A/pix (code 2); 1800 ln/mm 
grating, range 6200-6900 A, 2x2 binning, and 0.64 A/pix 
(code 4)'i an d 1800 ln/mm grating, range 4550-5200 A, 
2x2 binning, and 0.83 A/pix (code 10). 

Low and medium resolution, absolutely fluxed spec- 
tra of NCyg06 were also obtained with the AFOSC im- 
ager+spectrograph mounted on the 1.82m telescope oper- 
ated in Asiago by the INAF Astronomical Observatory of 
Padova. It is equipped with a Tektronix TK1024 thinned 
CCD, 1024x1024 array, 24 ^im pixel, with a scale perpen- 
dicular to dispersion of 0.67 "/pix. All observations were 
obtained with a 1.26" slit always aligned east-west. Four 
different instrument set-ups were used: 300 ln/mm grism, 
range 3500-7780 A, lxl binning, and 4.24 A/pix (code 6 
in Table 2); 1720 ln/mm volume phase holographic grism, 
range 6400-7050 A, lxl binning, and 0.64 A/pix (code 8); 
1280 ln/mm volume phase holographic grism, range 8270- 
9270 A, 1 x 1 binning, and 0.89 A/pix (code 9); and and the 
combination of a 79 ln/mm Echelle grism +150 ln/mm 
cross-disperser grism covering the whole wavelength range 
3600-8800 A in 13 orders with no gaps, providing a 3 000 
resolving power with lxl binning and 2.4 " slit width 
(code 7). 

High resolution spectra were obtained with the Echelle 
spectrograph mounted on the 1.82m Asiago telescope 
(code 5 in Table 2). It is equipped with an EEV CCD47-10 
CCD, 1024x1024 array, 13 /im pixel, covering the interval 
AA 3600—7300 A in 32 orders (wavelength gaps between 
adjacent orders are present redward of 4900 A). The slit 
was always aligned east-west, and its 2" projected width 
provided a 20 000 resolving power. 
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Fig. 2. Left panel: a small portion of the combined Echelle 
spectra for April 13 and 14, 2006 of Nova Cyg 2006, cen- 
tered on the interstellar Nal Di,2 absorption lines. Right 
panel: Blowup of the interstellar Nal D2 line (thicker line) 
showing the fit with two components at —39.1 (±0.4) and 
—11.1 (±0.3) km sec -1 (thinner line). 

3. Astrometry and progenitor 

Astrometry was performed on SRO images using SLALIB 
(Wallace 1994) linear plate transformation routines in 
conjunction with the UCAC2 reference catalog. Errors 
in coordinates were typically under 0.1 arcsec in both 
coordinates, referring to the mean coordinate zero point 
of the reference stars in each field. The coordinates we 
derived for NCyg06 are 

012000 = 21 /l ll m 32 s .33 (±0".12) 
5.72000 = +44°48'03".74 (±0".03), 

close to the coordinates measured by Yamaoka (2006b) at 
32 s . 35 (±0".15) and 03".66 (±0".14). The corresponding 
galactic coordinates are 1=087°. 372 and b=— 02°. 357. 

At this position, Steeghs et al. (2006) identified a faint 
point source on IPHAS images for Aug 3, 2004, measured 
at r'=20.30 (±0.05) and i'=19.76 (±0.07) and with ap- 
preciable Ha emission, which they argue is the progenitor 
of the nova. A search of archival Asiago Schmidt plates for 
1962-1984 was carried out by Jurdana-Scpic and Munari 
(2006) who found no brightening above the i?~18.5 plate 
limit. Both POSS-I and -II images also do not show a 
counterpart at the position of the nova. 

3.1. New Asiago archive plates 

In addition to the photographic plates already measured 
by Jurdana-Sepic and Munari (2006), we identified in the 
plate archives of the Asiago Schmidt telescopes another 
block of 57 direct-image plates covering the position of 
the nova progenitor over the period 1989-1997. All them 
have been inspected to search for a possible brightening 
during quiescence that could have caused the progenitor 
to rise above the i3~18.5 plate limit. None has been found. 
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Fig. 3. Comparison between the equivalent width of three 
diffuse interstellar bands and interstellar reddening for 
field stars measured by Megier et al. (2005). The crosses 
represent the position of Nova Cyg 2006. 

Table 3 (available electronic only) reports details on these 
plates. 

To estimate the .B-band magnitude expected in qui- 
escence for the nova progenitor, we adopted the mean 
fluxed spectrum from Zwitter and Munari (1995) for a 
cataclysmic variable in quiescence, increased its redden- 
ing to match the i?s_\/=0.56 affecting NCyg06 (see Sect. 
4), and scaled its flux to match the Steeghs et al. (2006) r' 
and i' IPHAS magnitudes. This caclulation indicates that 
the progenitor should be B~21.0 mag in quiescence. The 
negative detections on Asiago archive plates therefore ex- 
clude CV-type outbursts with an amplitude greater than 
A£?^2.5 mag. Frequent and large amplitude outburst ac- 
tivity as displayed by SS Cyg (<Ay>«4 mag, outbursts 
every «55 day, of «15 day duration; cf. Cannizzo and 
Mattei 1992), a prototype cataclysmic variable, should 
have been easily detected on the Asiago archive plates. 

3.2. The faint optical companion 

Seven arcsec west of the nova lies a faint field star, identi- 
fied with 2MASS 21113188+4448032 (J=13.716 ± 0.029, 
i?=13.441 ± 0.046, if s =13.480 ± 0.037). We measured 
on NOFS images for JD=2454247 its optical magnitude 
and colors as V=15.22, U - 5=+0.33, B - V=+0.70, 
V-R C =+QA2, i? c -/c=+0.46, V-J c =+0.88, with errors 
not exceeding 0.02 mag. Its astrometric position (J2000) 
from these images, with respect to the UCAC2 reference 
catalog, is a=21 11 31.88 (± 0.1 arcsec), <5=+44 48 03.2 
(± 0.1 arcsec). Steeghs et al. (2006) measured r'=14.8 for 
this star, in line with our i?c=14-80. 

The faint optical companion could not be resolved un- 
til the advanced decline set in ^260 days past maximum, 
when its contribution to the measurement of NCyg06 was 
no longer negligible. At later dates, the companion was 
spatially resolved only by instruments a and g thanks to 
better seeing and sampling. Photometry in Table 1 ob- 
tained with instruments b, c, d, e, and / always includes 
the contribution by the faint optical companion, which 
can be easily accounted for by using the companion's 
U BV Role magnitudes. The photometry in Table 1 ob- 
tained with instruments a and g only includes the contri- 
bution by the faint optical companion for t<260. At later 



dates, its was possible to measure the companion and the 
nova via PSF-fitting separately, and the values reported in 
Table 1 for i>260 and instruments a and g pertain to the 
nova alone. All values quoted in the text and the light- 
and color-curves presented in this paper (Figs. 1, 4, 10, 
15-20) have the contribution by the companion removed. 

4. Reddening 

Sharp Nal D12 interstellar absorption lines arc prominent 
on the high-resolution Echelle spectra of NCyg06 obtained 
on April 13 and 14, 2006, close to the optical maximum. 
Figure 2 shows the overall appearance of the interstellar 
lines superimposed onto the emission component of the 
wide Nal P-Cyg profile of the nova itself, as well as an 
enlargement of their profile that reveals how the interstel- 
lar lines are actually a blend of two distinct components 
of similar intensity. Their heliocentric radial velocities are 
-39.1 (±0.4) and -11.1 (±0.3) km sec" 1 , and their equiv- 
alent widths (for the D2 component at 5890 A) are 0.484 
and 0.487 A, respectively. Following the relation between 
the interstellar reddening and equivalent width of Nal D2 
line calibrated by Munari and Zwitter (1997), both com- 
ponents correspond to a reddening of 0.28 mag, for a total 
reddening of E\%= 0.28 + 0.28 = 0.56, which will be 
adopted throughout this paper. 

van den Bergh and Younger (1987) derived a mean 
intrinsic color (B — V) =— 0.02 ±0.04 for novae at t^. 
Figure 1 and the photometry in Table 1 show that NCyg06 
at ti had B — 1^=0.54. This corresponds to i?B_y =0.56, a 
perfect match to the result from interstellar lines. Russell 
et al. (2006) preliminarily estimated i?B_y=0.59 from in- 
frared spectroscopy. 

There are several diffuse interstellar bands (DIBs) vis- 
ible on the NCyg06 spectra. On the April 8 spectra (when 
the continuum was strongest and smoothest), DIBs at 
6203, 6270, and 6376-79 A are present but too weak 
for a meaningful measure against the overlapping P-Cyg 
components of Fell multiplet 74 lines. DIBs at 4430, 
5780, 6284, and 6614 A show equivalent widths of 2.57, 
0.75, 1.17, and 0.14 A, respectively. A comparison of the 
equivalent width of the last three DIBs with the redden- 
ing/equivalent widths measured for field stars by Megier et 
al. (2005) is presented in Fig. 3. The correlation between 
reddening and DIB equivalent widths is known to be a 
loose one (cf Herbig 1995). Nevertheless, it is worth notic- 
ing that, while the intensity of DIB 6614 A in NCyg06 
seems to nicely follow the mean relation for field stars, 
the DIBs at 5780 and 6284 are appreciably stronger than 
in field stars. DIB 4430 A similarly lies above the mean 
relation for field stars (though still within the scatter of 
individual data) presented by Isobe et al. (1986). Overall, 
the DIBs toward NCyg06 are somewhat stronger than ex- 
pected for the Eb-v=0-56 reddening affecting the nova. 
The nova itself should not be contributing to the observed 
DIBs because their very low observed heliocentric radial 
velocities are what is expected for an origin in the inter- 
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stellar medium and bear no relationship to the high radial 
velocity of the ejecta. 

5. Early evolution 

The early evolution of NCyg06 was that of a fairly normal 
Fell-type nova, from both the spectroscopic and photo- 
metric points of view, with no hint of the later peculiari- 
ties that led to one of the most peculiar lightcurves ever 
recorded for a nova. 

5.1. Photometry 

The first two months of the nova light-curve are presented 
in Fig. 4. The optical photometric evolution of NCyg06 
was characterized by a fast rise to a maximum of V=8.0 
(±0.1) on Apr 5.6 (±0.5) 2006 UT (=JD 2453831.1). 
The time required to raise the final two magnitudes was 
£ r! 2^2.3 days, and the nova was discovered 2.8 days be- 
fore maximum, when it was 2.5 mag fainter. A week before 
maximum, the nova was below the patrol plate limit at 12 
mag (cf IAUC 8698). The decline in V-band during the 
first two months (4<Ai<56, where At=t-t max is the time 
in days from maximum) strictly followed the exponential 
slope given by 

V(t) = 7.6 + 2.4 x log(Ai) (1) 

shown in Fig. 4 as a solid curve. A large (Af?«0.3 mag) 
variability of NCyg06, but with constant color (AB — 
V~0.0), is superimposed on the general decline during the 
first ten days after maximum. This scatter cannot be as- 
cribed to observational problems. In fact, the correspond- 
ing observations have the highest accuracy, and the results 
are confirmed from independent instruments, as well as re- 
peated observations with the same instrument. 

The nova decline times were tjf =10.4 an d t^f =24 
days, on the borderline between very fast and fast novae. 
Both t2 and £3 are in normal proportion among them. 
In fact, the following relation fits nicely the sample of 
~20 best-studied novae away from the galactic Bulge (i.e. 
45°<galactic longitude<305°), 

h ~ 6 x t° 2 7 - 7, (2) 

and NCyg06 falls right on it. The rise t T< 2 and decline £2, ^3 
times of NCyg06 are close to the proportionality relation 
found by Schmidt (1957): log t T , 2 = -0.3 + 0.7 log t 2 
= —0.5 + 0.7 log £3. In addition, the ~13.3 mag outburst 
amplitude and i5f =10.4 fall precisely on the Warner (1995) 
relation between amplitude and t 2 . 

van den Bergh and Younger (1987) found that the aver- 
age intrinsic color of novae at maximum is (B— V^) o =+0.23 
±0.06 (dispersion 0.16 mag). Our first color measurement 
was obtained close to maximum at At~+l d .6 (cf. Table 1). 
It gave B — V=+0.85, which corrected for the Eb-v =0.56 
reddening, corresponds to (B — V) o =+0.29, thus reason- 
ably close to the mean of intrinsic colors of novae at max- 
imum. 




HJD since discovery on 245382B.307 



Fig. 4. The early photometric decline of Nova Cyg 2006. 
For instrument identification see Sect 2.1. The solid curve 
over-plotted to the V data in panel d is Eq. (1), while 
the dotted curve is drawn by hand only for eye guidance. 
The a panel shows the evolution in equivalent width of an 
un-blended emission line for each of the Fell multiplets 
27, 42, 49, 55 and 74 (4233, 4924, 5317, 5535 and 6239 A, 
respectively). Panel b presents the change in radial velocity 
of absorption systems (the "de" diffuse-enhanced and the 
"pi" and "^2" principal systems seen in Ha, and the "Or" 
Orion system seen in Hel 6678 A). 

5.2. Distance and bolometric magnitude 

Both the rate of decline and the observed magnitude 15 
days past maximum are popular methods of estimating 
the distance to the nova. 

Published relations between absolute magnitude and 
the rate of decline generally take the form M max = 
a n log t n + f3 n . Given the lack of observations in B band 
right at maximum light, we are limited to equations valid 



U. Munari et al: Nova Cyg 2006 7 



At 



36.4 



21.5 



18.4 



6.5 



J2.5 




4000 



5000 



6000 



7000 



Fig. 5. Early spectroscopic evolution of Nova Cyg 2006. The spectra are continuum-normalized and shifted to avoid 
overlap. The ordinates are in logarithm of the flux to emphasize the visibility of weaker features. At and AV are 
days and magnitudes from maximum. Emission lines are identified (Fell multiplets by their numbers and comb-like 
markings). 



for V band alone. Using the Cohen (1988) V-ti relation, 
the distance to NCyg06 is 8.0 kpc, and it is 7.3 kpc accord- 
ing to the Schmidt (1957) V-t 3 relation. Using the dclla 
Valle and Livio (1995) expression calibrated on novae ap- 
pearing in the LMC and M31, the distance to NCyg06 
would increase to 9.7 kpc. Since it is not applicable to 
galactic novae, it will be not considered further. 

Buscombe and de Vaucouleurs (1955) suggested 
that all novae have the same absolute magnitude 15 
days after maximum light. There are several calibra- 
tions of it: MY 5 =— 5.2±0.1 from Buscombe and de 
Vaucouleurs (1955), M^=-5.60±0.43 from Cohen 
(1985), M^=-5.23±0.16 from van den Bcrgh and 
Younger (1987), M^=-5.38 from van den Bergh 

(1988) , M^=-5.69±0.14 from Capaccioli et al. 

(1989) , Mff=-5.85 according to Schmidt (1957), 
Mff=-5.50±0.18 from de Vaucouleurs (1978), and 
Mf 5 =-5.74±0.60 from Pfau (1976). The brightness of 
NCyg06 15 days after maximum light was Ui 5 = 10.45, 
.Bi5=10.97. Comparing with the above absolute mag- 
nitudes and correcting for i?s_y=0.56, the distance 
turns out to be 6.1, 7.3, 6.2, 6.6, 7.6, 8.0, 6.8, and 7.6, 
respectively (approximating m pg KiB). 

Taking the unweighted mean of the ten determinations 
above, the distance to NCyg06 is d—7.2 kpc, with an error 
of the mean of 0.2 kpc and a dispersion ad = 0.65 kpc. The 
bolometric correction for an expanded photosphere with 
B - U=+0.29 is B.C.=-0.02 (from Drilling and Landolt 



2000). Coupled with a 7.2 kpc distance, i?B_y=0.56 and 
U max =8.0 mag, it gives for NCyg06 a bolometric magni- 
tude at maximum of Mboi=— 8.0 mag. 

At NCyg06 galactic latitude, a 7.2 kpc distance trans- 
lates to a height above the galactic plane of z=0.3 kpc, still 
within the vertical scale height of the galactic Thin Disk. 
At galactic longitude ^=87°. 37, the line of sight to NCyg06 
crosses the Perseus spiral arm 5 kpc away from the Sun, 
and stops midway between the Perseus and Cygnus spi- 
ral arms (cf. Vallee 2005). According to the Brand and 
Blitz (1993) maps, the radial velocity of interstellar mat- 
ter belonging to the Perseus spiral arm where it is crossed 
by the line of sight to NCyg06 is w— 40 km sec -1 , so the 
—39 km sec" 1 component seen in the high-resolution pro- 
files of interstellar Nal D12 (Fig. 2 and Sect. 4) is quite 
possibly associated with the Perseus spiral arm. The other 
component at —11 km sec -1 should be connected to ma- 
terial in front of the Perseus spiral arm, as suggested by 
the Neckel and Klare (1980) extinction maps and by the 
lower observed radial velocity. 



5.3. Spectroscopy 

The low-resolution, broad-wavelength range spectroscopic 
evolution of NCyg06 during the first two months of the 
outburst is presented in Fig. 5. 
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Fig. 6. Evolution of the Ha profile at three epochs during 
the first three weeks past maximum (see Sect 5.3.1 for 
details). 



5.3.1. Absorption systems 

The evolution of the Ha profile during the first three 
weeks of the outburst is shown in Fig. 6. It is based on 
high-resolution observations obtained at epochs At=+2.5, 
+6.5, and +18.4 days. A fourth profile for At=+21.5, sim- 
ilar to that of day +18.4, is not shown to save space, but 
it is nevertheless included in Table 2, which reports the 
results of the Ha profile deconvolution at all four epochs. 



Table 4. Heliocentric velocity, velocity span at half max- 
imum, equivalent width, and integrated absolute flux (in 
units of 10~ 13 erg cm' 2 sec' 1 ) of the emission component 
and pre-maximum, principal and diffuse enhanced absorp- 
tion systems for the Ha profiles shown in Fig. 6 (and in 
addition for the Ha profile of Apr 27.09, 2006 not shown 
in Fig. 6 to save space). 
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We do not consider later observations because of the re- 
inforcement of [Nil] emission in May and the blending of 
its 6548, 6584 A lines with the Ha. 

The Ha profile fitting was carried out assuming that 
the underlying emission component has a symmetric pro- 
file and that absorptions have Gaussian profiles. The pro- 
file of the emission component is therefore derived from 
the portion of the observed profile that extends from the 
peak to the asymptotic merging of the red Ha wing into 
the continuum. In Fig. 6 the profile of the emission compo- 
nent is plotted with a dashed line, those of the absorption 
components by solid lines. 

The first spectrum (At— +2.5), obtained when the 
nova had already declined by AV=0.7 mag, shows four 
absorption components, listed in Table 4. They corre- 
spond to the simultaneous presence of the classical pre- 
maximum, principal, and diffuse- enhanced absorption sys- 
tems (indicated by pm, p, and de in Fig. 6, respectively), 
which are common to most novae and are described in 
detail by McLaughlin (1960, hereafter McL60). A note- 
worthy characteristic of NCyg06 is the splitting of the 
principal absorption spectrum into two components (pi, 
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p 2 ), as seen previously in other novae (e.g. N Aql 1918 
and N Gem 1912). At AV—0.7 mag, the pre-maximum ab- 
sorption system was already declining in equivalent width, 
while the principal system was rapidly gaining strength 
and the diffuse- enhanced system had just emerged. Their 
heliocentric radial velocity, velocity width at half intensity 
and integrated absolute flux are given in Table 4. 

Data summarized by McL60 show a correlation be- 
tween mean radial velocity of the various absorption sys- 
tems and the speed class of the nova (mean velocity be- 
cause the velocity of any given system usually changes 
with time). This correlation states Wpro-max^— 4750/i 2 for 
the pre-maximum absorption system. The t 2 —lOA days 
for NCyg06 would correspond to —450 km/s, close to the 
observed —400 km/s (Table 4). The velocity relation for 
the principal system is logw pr i n = 3.70 — 0.51ogi 3 = 3.57 

— 0.51ogt 2 , and predicts w-1080 km/s for NCyg06. The 
epoch means of the pi,P2 velocities in Table 4 (weighted 
according to their equivalent widths) are —850, —1087, 

— 1103, and —1112 km sec" 1 . Their global average is 

— 1040 km sec -1 , again close to predictions. The McL60 
relation for the diffuse enhanced system is logwdif-cnh 
= 3.81 - 0.41ogt 3 = 3.71 - 0.41ogi 2 , and it predicts 
«— 1880 km/s velocity for NCyg06, in good agreement 
with the averaged observed value of —1830 km/s. 

All three absorption systems in NCyg06 appear to 
move ~50 km sec -1 more slowly then predicted on the 
basis of the relations by McL60. The difference is actually 
more than twice as large. In fact, the component of the 
radial velocity due to the galactic rotation can be written 
as 

RVq — w sinfe + u cosl cosh — v sin/ cos6 

- 2[A(R - R ) + a(R - R ) 2 } sin? cosfe (3) 

where R 2 = R 2 +d 2 - 2R dcosl, A = -(R /2)(du/dR) Ro 
and a = -(R /4:)(d 2 uj/dR 2 ) Ro . Here, R and R a are 
the nova and Sun galacto-centric distances, d is the 
nova-Sun distance, (l,b) the nova galactic coordinates 
and (u 0} v ,w ) is the solar motion vector. Adopting 
(u 07 v ,w ) = (10.0,5.3,7.2) km sec -1 from Dehnen and 
Binney (1998), and A = 17.0 ± 1.5 km sec -1 kpc -1 and 
a = -2.0 ± 0.6 km sec -1 kpc -2 from Hron (1987), the 
radial velocity of NCyg06 expected from galactic rotation 
would be —69 km sec -1 . Therefore, it seems that the mean 
barycentric outflow velocity of the pre-maximum, prin- 
cipal and diffuse-enhanced spectra of NCyg06 are about 
120 km sec -1 slower than predicted by the McL60 rela- 
tions. 

The Fell line profiles in the low-resolution spectra for 
April 8 and 12 in Fig. 5, and the Echcllc spectra for April 
13 and 14, also show the presence of the principal and 
diffuse-enhanced absorption systems. The pre-maximum 
system had become too weak to be detected on the Fell 
lines by the time of the Echelle spectra. The radial ve- 
locities of the absorption systems on the Fell lines were 
close to what was measured for Ha. On the low-resolution 
April 8 spectrum, they were ~— 830 for p\ principal and 



Table 5. Radial velocity and equivalent width of Orion 
the absorption system as observed in the Hel 6678 A line. 
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56.5 


3.8 
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w— 1700 km/s for diffuse-enhanced. On the Echcllc spec- 
trum for April 12 they were «— 985 for pi principal, 
rs— 1295 for pi principal, and w— 1780 km/s for diffuse- 
enhanced. Also the ratio of the equivalent widths of p\ 
and diffuse-enhanced was similar in Ha and Fell lines: 
1.29 and 1.30, respectively, on April 12. The only major 
difference with Ha was in the much weaker intensity in the 
Fell lines of the p 2 system compared with the p\ system. 

The time behavior of the radial velocity of p\ , p 2 prin- 
cipal and diffuse-enhanced absorptions on the Ha profile 
is illustrated in Fig. 4. The velocity of the p\ system has 
been increasing with time as t>ioc— 215 log(Ai), and for p 2 
as u 2 oc— 4001og(At) km/s. The velocity increase was min- 
imal for the pre-maximum absorptions, while there has 
been no increase at all in the diffuse-enhanced systems. 
This agrees with what is seen in most novae (cf. McL60). 

Feeble traces of the emergence of the Orion absorp- 
tion system appear on the April 12 spectrum (Ai=+6.5, 
AV=1.6 mag) of Fig. 5. The Orion absorptions grew 
much stronger on the April 24 and 27 spectra (Ai=+20, 
AV=2.8 mag). It was still weakly present on the May 11 
spectrum (At=+36, AV=3.4 mag). The Orion system af- 
fected the Hel, Nil and Oil lines but not the hydrogen 
ones, as is common in many novae. The radial velocity 
and equivalent width of the Orion absorption component 
of Hel 6678 A are listed in Table 5, while Fig. 4 shows 
its velocity progression in comparison with the other ab- 
sorption systems. In many novae, the velocity of the Orion 
absorption is comparable to or higher than for the diffuse- 
enhanced system, while the Orion system moves about 
250 km/s slower in NCyg06. 

In typical novae, as the numerous Nil and Oil lines 
of the Orion system fade away, a pair of NIII absorption 
lines emerges at 4097 and 4103 A. This also occurred in 
NCyg06, for which the NIII absorption is quite visible on 
the blue wing of the HJ emission line of the June 1 spec- 
trum (At=+56, AV=3.8 mag) illustrated in Fig. 9. Its 
velocity is ~— 2200 km sec -1 . The June 1 spectrum cor- 
responds to the date when the nova abandoned the ex- 
ponential decline of Eq.(l), began to slow its decline rate 
toward the plateau phase of mid 2006 and then rose to 
the second maximum of late 2006. This anomalous photo- 
metric behavior was mirrored by an initially slowed, and 
then aborted, spectroscopic evolution. In fact, by the time 
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the NIII 4100 A absorption appears, a typical nova is en- 
tering the transition phase, when its spectrum turns from 
the stellar to the nebular type. This did not happen in 
NCyg06. It did not go through the transition phase and 
the development of the NIII absorption soon stopped. The 
equivalent width of the latter, in fact, reached a modest 
1.1 A on the June 1 spectrum and did not rise to the 
expected ^10x larger equivalent width it should have at- 
tained at later dates. Instead, it rapidly faded away, leav- 
ing no obvious trace on the June 22 spectrum in Fig. 9. 

5.3.2. Emission lines 

Figure 5 shows the evolution of the emission lines dur- 
ing the first two months of the outburst. On At=+2.5, 
the emission lines were essentially restricted to the hy- 
drogen Balmer series and emerging Fell lines (principally 
from the multiplets 42, 49, and 74), on top of a poorly 
structured, relatively cool continuum. The intrinsic color 
at that time, (B-V) o =+0.29, was still only marginally per- 
turbed by emission lines and corresponded to that of an 
A9/F0 III star. The emission line spectrum was typical of 
the Fe77-class defined by Williams (1992). The expansion 
velocity (corrected for instrumental resolution) derived 
from the width at half maximum was similar for all hydro- 
gen Balmer emission lines and amounted to 780 km sec -1 ; 
the velocity estimated from the width of Fell lines is 
25% slower. Four days later, on At=+6.5 the emission 
line spectrum was getting much stronger (cf. the increase 
in equivalent width of Ha and Fell lines in Figs. 1 and 
4), and the Fe//-classification even clearer. Auroral [Nil] 
5755 A and [01] 5577 A emission lines were becoming vis- 
ible, as well as nebular [OI] 6300, 6364 A. The intensity 
of the auroral lines was comparable to that of the nebular 
transition lines, indicating high densities in the outermost 
neutral external regions of the ejecta where they formed. 

On Ai=+18.4, AV=2.6 mag, the [OI] lines grew in in- 
tensity more than the other lines, reminiscent of the [OI] 
flash (cf. PG57) seen in other novae at a similar evolu- 
tion stage. The same happened to the [Nil] lines and the 
NII+OII+NUI blend at 4640 A on Ai=+36.4, AV=3.3 
mag when Fell weakened and Hel turned into weak emis- 
sion, corresponding to the [Nil] flash and 4640-phase ob- 
served in many novae at A 1^=3. 0/3. 5 mag (cf. PG57). 

The reduced intensity of auroral with respect to nebu- 
lar [OI] transitions in May indicates a decline in electronic 
density in the outer ejecta, as expected by the ongoing ex- 
pansion and persistently low ionization stage. 

5.3.3. 01 lines 

Emission lines from neutral oxygen offer several inter- 
esting applications for novae, some of which have been 
discussed by Strittmatter et al. (1977), Williams (1994), 
Bhatia & Kastncr (1995), Kastner & Bhatia (1995), and 
Osterbrock (1989). Formalism and atomic constants in the 
present section were adopted from these sources. 



Nova Cyg 2006 

The [OI] 5577, 6300, 6364 A lines are among the first 
forbidden lines to appear in nova spectra thanks to their 
high critical densities. They persist even when [OIII] dom- 
inates the emission spectrum and high ionization lines ap- 
pear (see for example the simultaneous presence of [FeVII] 
and [OI] lines in the spectra of Fig. 14). In the simpli- 
fied case of spherically symmetric, homogeneous ejecta, 
created by a steady wind and photoionizcd by a central 
source, the [OI] lines trace ionization-bounded conditions 
and are expected to form in the external, neutral regions. 
The ratio of the 6300, 6364 A nebular transitions under 
optically thin conditions is 3:1 from their transition prob- 
abilities, and it scales with optical depth r (in the 6300 A 
line) as 



At typical oxygen abundance, ejected mass and outflow 
velocity for a nova, [OI] 6300 A should turn optically 
thin within the first few days after maximum brightness. 
Williams (1994) noted how the 6300/6364 A flux ratio in 
novae is almost always lower than 3:1, indicating persis- 
tently optically thick conditions in the lines. He proposed 
that optically thick [OI] lines in novae come from small, 
very dense, neutral globules embedded in ambient, ionized 
ejecta and that the globules are internally powered by the 
decay of unstable isotopes produced by initial nuclear re- 
actions (e.g. Hernanz 2005, and reference therein, for a 
discussion of unstable nuclides). 

Figure 1 illustrates the evolution of 6300/6364 A flux 
ratio over the whole outburst of NCyg06. The lines first 
appeared on the April 12 spectrum, ~10 days after the 
outburst onset, with a 2.2:1 ratio, and then initially 
evolved toward the 3:1 optically thin value. However, by 
late May 2006, when the decline of the nova slowed, the 
evolution of the 6300/6364 A flux ratio reversed its course 
and progressively decreased, reaching the minimum value 
1.1:1 by the time of second maximum. Such a low value 
requires r «7. As soon as the nova entered the transi- 
tion phase immediately following the second maximum, 
with the ionization spreading throughout the ejecta, the 
6300/6364 A ratio rapidly converged to the 3:1 value ex- 
pected for optically thin conditions. The 6300 A integrated 
line flux, on the other hand, displayed only a modest de- 
cline during the transition. Such a decline was confined to 
the phase when optically thin dust condensed (cf. Sect. 7 
below), and afterwards the 6300 A line progressively and 
rapidly regained the flux level it had reached before the 
transition phase. 

The [OI] lines provide an estimate of the mass in 
the neutral region of the ejecta. Table 6 reports the re- 
sults for some distinct epochs in the evolution of NCyg06, 
where Eq. (4) is used to derive the optical depth from 
the 6300/6364 A flux ratio. The high-density conditions 
in the neutral regions make the auroral 5577 A line easily 
detectable and therefore usable with the nebular 6300 A 
line to estimate the electron temperature. How the lat- 
ter correlates with the (reddening corrected) 5577/6300 A 
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Table 6. The optical depth, electron temperature and 
mass of neutral oxygen in NCyg06 from the [01] 5577, 
6300 and 6364 A emission lines. The listed epochs are 
those of the spectra showing the 5577 A line in emission. 



At 

(days) 



-^6300 
-F6364 



76300 



-^6300 
-F5577 



^"elec 

(°K) 



Moi 
(M Q ) 



6.5 


1.304 


4.23 


2.133 


4316 


1.2e-ll 


8.1e-05 


18.4 


2.290 


0.90 


1.189 


5927 


2.3e-ll 


1.3e-05 


21.5 


2.479 


0.61 


1.330 


5948 


2.2e-ll 


1.0e-05 


36.4 


2.474 


0.62 


3.498 


4859 


1.4e-ll 


1.7e-05 


50.3 


2.417 


0.70 


4.336 


4640 


1.2e-ll 


1.7e-05 


56.5 


2.288 


0.90 


5.857 


4338 


1.5e-ll 


3.5e-05 


78.3 


2.190 


1.07 


10.844 


3888 


1.0e-ll 


4.5e-05 


79.5 


2.359 


0.79 


11.503 


3924 


1.2e-ll 


5.0e-05 


96.3 


2.283 


0.91 


10.412 


3952 


1.0e-ll 


3.9e-05 


270.1 


2.548 


0.52 


4.764 


4629 


4.4e-12 


6.1e-06 


282.2 


2.553 


0.51 


4.243 


4731 


4.9e-12 


6.1e-06 


290.2 


2.570 


0.49 


4.108 


4768 


7.4e-12 


8.8e-06 


339.6 


2.796 


0.22 


6.138 


4537 


5.0e-12 


6.7e-06 



At: time since V brightness (April 5.6, 2006) 
t"63oo: optical depth in the 6300 A line 
T e i ec : electronic temperature 
Moi: mass of neutral oxygen 

The flux ratios (F e300 /F e364 and ^6300/^5577) and the 
^6300 integrated line flux (in erg cm~ 2 sec -1 ) are 
corrected for the Eb-v=0.56 reddening affecting NCyg06 



flux ratio for a wide range of 6300 A optical line depths 
is shown in Fig. 7. The proportionality between the inte- 
grated flux in the 6300 A line and the amount of neutral 
oxygen in the ejecta as a function of optical depth and 
electron temperature is highlighted in Fig. 8 for the 7.2 
kpc distance and i?e_y=0.56 reddening derived earlier for 
NCyg06. The high density condition ruling the region of 
formation of OI lines is highlighted by the OI i*8446 / i^6300 
emission line ratio. Its reddening-corrected value was 16.5 
on Ai=18.4 and Ai=21.5. Comparing with the computa- 
tion by Kastner and Bhatia (1995), it suggests a density 
> 10 10 cm~ 3 . The ratio lowered to 0.4, 2.8, 0.4, and 0.2 
on Ai=186, 253, 340, and 527, respectively, indicating a 
density > 10 8 cm~ 3 . 

The data in Table 6 suggest that, at the beginning of 
the nova evolution, the neutral oxygen was present pre- 
dominantly in the diffuse ejecta, while at the later stages 
it only survived in isolated, dense, and neutral globules. 
The mass of these globules accounts for <~ 1/12 of the total 
mass of the ejecta. This rough proportion can be derived 
by comparing the mass at first [OI] detection (At=6.5 
days), when nova ejecta were predominantly neutral, with 
the mass at advanced decline (Ai>260 days) when the 
nova ejecta were essentially completely ionized (density 
bounded). 



It will be shown in Sect. 11 that the ejecta of NCyg06 
are oxygen enhanced (an O/H ratio ^32 x larger than 
in the Sun). On Ai=6.5 days, all oxygen in the NCyg06 
ejecta was neutral, and therefore the mass derived for neu- 
tral oxygen corresponds to the total mass of oxygen in the 
ejecta (Mo=Moi). Noting that the oxygen mass fraction 
derived in Table 9 from nebular analysis is 0.163 and that 
dust had not yet formed during early decline, the total 
mass of the ejecta can be estimated as 



M, 



ejecta — 



0.163 



5 x 10" 



Mr. 



(5) 



where M O i=8.1xl0~ 5 M Q is taken from Table 6. An in- 
dependent estimate for M e j ecta will be obtained in Sect. 11 
from analysis of the emission line spectrum in the nebular 
phase. 

It is also noteworthy in Table 6 that the amount of neu- 
tral oxygen parallels the photometric evolution. In fact, 
in Sect. 5.1 we discussed how by June 1, 2006 the nova 
was departing from a pure exponential decline (to pro- 
gressively slow its decline and later rise to a second max- 
imum). At the same time, the amount of neutral oxygen 
(and presumably of neutral ejecta) was increasing, as if 
the ionization boundary was retracting (in a mass sense) 
through the ejecta. 

The intensity of the OI 8446 A emission line under 
normal recombination, optically thin conditions should be 
3/5 of the OI 7774 line. The reddening-corrected flux ratio 
observed in NCyg06 is F 8446 / F 7774 =7 .7 on Ai=18.4 (de- 
cline from first maximum), 21.5 on Ai=290 (decline from 
second maximum), and 12.5 on At=340 (nebular condi- 
tion during the advanced stages of the final decline). The 
inversion in intensity between the two OI lines is usually 
associated with fluorescence pumped by absorption of hy- 
drogen Lyman-/3 photons, as first pointed out by Bowcn 
(1947). For the Lyman-/3 fluorescence to be effective, the 
optical depth in Ha should be large, presumably owing 
to the population of the n = 2 level by trapped Lyman-a 
photons. The ^8446/ Fua under optically thin, low ioniza- 
tion conditions and typical nova chemical abundances is 
quite low, ~10~ 3 (Strittmatter et al. 1977). The evolu- 
tion of the i ? 8 44 6 /Ff/ Q flux ratio for NCyg06 is presented 
in Fig. 1. It reached a much higher ^0.5 value during 
early evolution (on Ai=18.4 and 21.5). It indicates a very 
large optical depth in Ha, in agreement with the ejecta 
still being predominantly neutral and very dense at that 
time. At later times the picture was more confused due 
to the dumpiness of Ol-cmitting regions described ear- 
lier and the growing [Nil] contribution to the Ha ob- 
served flux. Nevertheless, even at advanced evolutionary 
stages, the F^q/Fho did not decline below 0.05, indi- 
cating persistent (even if reduced) optically thick condi- 
tions in Halpha. This agrees with the evidence from the 
Fs446 / -F7774 ratio that Lyman- (3 fluorescence has remained 
effective throughout all NCyg06 evolutionary stages. 
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Fig. 9. Spectroscopic evolution of Nova Cyg 2006 during the plateau phase in June-September 2006 (see Fig. 5 for 
plotting details). The asterisks mark the positions of the principal Orion absorption lines on the September 26 spectrum 
(see Sect. 6 for the line list). 



6. The plateau phase 

When NCyg06 left the exponential decline of Eq.(l) 
around June 1, it had magnitudes and colors of V=11.83, 
B - ^=+0.28, V - fl c =+1.44, R c - 7 c =+0.20 and 
V — Jc=+l-64. A minimum brightness was reached on 
July 15 at V=12.15, B - V=+0.15, V - i? c =+1.28, 
R c - J c =-0.01 and V - J c =+1.27. During the follow- 
ing rise in brightness, the same TZ-band brightness as for 
June 1 was regained 111 days later on September 20, 
when it was V=11.83, B - V=+0.16, V - i? c =+0.91, 
Rc - I C =+0.07 and V - 7 c =+0.99. 

The significantly bluer colors on September 20 com- 
pared to June 1 indicate that the nova had appreciably 
evolved. This is clear from the evolution of the spec- 
trum (illustrated in Fig. 9) during this plateau period. 
Progressively, Fell emission disappeared, while NII/H/3 
and NIII+H5/H7 increased, and [Nil] decreased with re- 
spect to Nil and increased compared to [OI]. However, 
the most interesting feature of this plateau period is the 
re-emergence of the Orion absorption system. The Orion 
absorptions were last seen on the May 11 spectrum and 
marginally still visible on June 1. On the August 21 and 
September 26 spectra of Fig. 9, the most notable Orion 
absorption lines are those of Hel 3889 A (—2140 km/s ra- 
dial velocity), 4471 (-2065 km/s), 5876 (-2260 km/s), 
Oil 4321 (-2125 km/s), Nil 4601 (-2100 km/s), 4788 
(-2150 km/s), 5672 (-2100 km/s). Particularly interest- 
ing is the re-emergence of NIII 4097-4103 (-1950 km/s) 



on the September 26 spectrum (compare the HS blue wing 
with that of the June 1 spectrum on Fig. 9). The average 
velocity of the Orion absorptions on September 26 was 
—2100 km/s, significantly higher than on June 1 when 
the Orion spectrum was last seen (cf. Table 5). The accel- 
eration seen in the Orion absorptions had the reverse sign 
for NIII 4097-4103: -2200 on June 1 and -1950 km/s on 
September 26. 



During this plateau period, we obtained high- 
resolution Ha profiles on June 1 and 24, July 18 and 
29, August 31 and September 26. The Ha profles from 
June 1 to August 31 are all the same and very similar 
to that of June 24 shown in Fig. 12, which is character- 
ized by a multi-peaked top and by 1340 km/s full width 
at half intensity. On the other hand, the September 26 
profile is characterized by the emergence of a hazy, broad 
and featureless pedestal whose full width at zero inten- 
sity approaches 6000 km/s, while the width at half inten- 
sity is ^4200 km/s (matching the ^—2100 km/s expan- 
sion velocity of the Orion absorption lines) . This pedestal 
grew conspicuously in relative intensity by the time of sec- 
ondary maximum (cf Fig. 12). It is worth noticing that the 
decline of the Orion absorptions in novae is usually accom- 
panied by the emergence in emission line profiles of hazy 
and broad components, even if Balmer lines do not always 
display them. 
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7. The second maximum 

The second maximum was reached by NCyg06 at V—9.95 
around HJD 2454071.4 (December 1, 2006) at At=+239.0 
and AV=1.95, with colors U - S=-0.29, B - ^=+0.28, 
V-R c =+0.68, i? c -^c=+0.28, y-/ c =+0.98. The emis- 
sion lines at that time had a minimal impact on the over- 
all energy distribution, and the intrinsic (B-V) =— 0.28 
corresponds to that of a B0 giant, compared to an F0 gi- 
ant at first maximum. Compared to the first maximum, 
the colors at the second maximum were bluer by the fol- 
lowing amounts: A(B - V)«-0.60, A(V - i?c)«-0.42, 
A{R C - 7c)«-0.20 and A(V - J c )«-0.52 (colors at t=0 
are extrapolated by parabolic fitting to the early decline 
data in Table 1 and Fig. 1). The time required to rise the 
last magnitude in the V band was 22 days, to decline the 
first two magnitudes was 9.6 days, and to decline the first 
three magnitudes was 11.9 days. The brightness and color 
evolution around the second maximum is highlighted in 
Fig. 10. 

The overall spectroscopic aspect of NCyg06 at the sec- 
ond maximum was appreciably different than at the first 
one, as seen when comparing the spectral evolution in 
Fig. 5 with Fig. 11). In addition to the Balmer series, 
only low-excitation absorption lines (mainly Call, Nal, 
Fell) were seen at the first maximum, in agreement with 
the F0 underlying spectral energy distribution. At second 
maximum, the underlying continuum was much hotter and 
the excitation of the absorption spectrum correspondingly 
much higher. The strongest, non-Balmcr absorption lines 
were: CII 6789, 6746, 5143; Nil 6482, 5678, 5667, 5045, 
5007, 4791, 4602; NIII 4639, 4513, 4100; Oil 4495, 4468, 
4416, 4318, 4070 and Hcl 6678, 5876, 4471, 3889. Their 
mean radial velocity was —2190 km/s for CII, —2135 km/s 
for Nil, -2070 km/s for Oil, and -2170 km/s for Hel. The 
absorption spectrum at second maximum represents the 
re-emergence of the Orion system that disappeared during 
the plateau phase. 

Quite interesting are the absorption components 
seen in Ha during the second maximum. The spec- 
trum for November 5, 2006 in Fig. 12 is represen- 
tative of them. The underlying emission is obviously 
much more structured compared to the simple, sym- 
metric one seen during first maximum (cf. Fig. 6), 
and therefore unambiguous identification of weaker ab- 
sorptions is severely affected. Nevertheless, at least two 
strong absorptions are safely identified: one centered at 
-905 km/s (VHM=145 km sec" 1 , e.w.=1.90 A, inte- 
grated flux=8.4x 10~ 13 erg cm" 2 sec -1 ), the other at 
-2215 km/s (VHM=355 km sec" 1 , e.w.=7.67 A, inte- 
grated flux=2.0xl0~ 12 erg cm" 2 sec" 1 ). They could be 
tentatively identified with the rc-cmergencc of princi- 
pal and diffuse- enhanced absorption systems seen during 
the decline from the first maximum (cf. Table 4). On 
November 22, the faster of the two components had grown 
to an integrated flux of 3.5xl0 -13 erg cm -2 sec -1 and an 
equivalent width 9.8 A. 



Similar to what has just been described for the ab- 
sorption spectrum, the emission spectrum at second max- 
imum was quite different from that at first maximum, 
reflecting the much higher temperature of the underly- 
ing continuum. At second maximum, Fell emission lines 
were not seen and were replaced by Nil, NIII, Oil, [OI], 
and Hcl. [Nil] 5755 disappeared, and the F 630 o /^6364 ratio 
approached 1.1:1 (cf. Fig. 1 and Sect. 5.3.3). Overall, the 
strong decrease in equivalent width of the emission lines at 
secondary maximum (cf. Ha behavior in Fig. 1) essentially 
stemmed from the increase in the underlying continuum. 
The evolution of integrated flux of emission lines went un- 
affected through the development of second maximum, as 
the evolution of [OI] and Ha in Fig. 1 illustrates. 

8. Final decline 

Soon after the second maximum was reached on December 
1, NCyg06 entered the final decline that is continuing up 
to now. During this phase, NCyg06 resumed the normal 
spectroscopic evolution of a nova that was interrupted by 
the plateau phase and rose to secondary maximum. The 
decline is characterized by a sharp change in the slope 
of the lightcurve occurring around December 12, as illus- 
trated by both Figs. 1 and 10. That date also marks the 
transition of NCyg06 from stellar to nebular spectral type. 

8.1. Decline prior to December 12 

Between the second maximum and December 12, the 
nova became fainter by AV—2.S mag in 11 days, with 
t^=9.6±0.2 and i^=11.9±0.3 days. The slope of the de- 
clines from the first and second maxima are quite different, 
as the respective £2/^3 suggest and a comparison of Figs. 4 
and 10 supports. The rate of decline from first maximum 
was slowing with time, that from second one was increas- 
ing with time. The B — V color initially turned bluer as the 
stellar continuum turned hotter, but this was short-lived 
and by December 8 it reversed the trend. By December 13 
it had resumed the smooth course it has been displaying 
since the earliest outburst phases (cf. Fig. 10). The V — Rc 
color increased almost exponentially from December 7 to 
12, as a result of the rapidly growing relative contribution 
of Ha to i?c _ band (cf. Ha equivalent width evolution in 
Fig. 1). The extra flux in the /c-band, clearly influencing 
Rc — Ic an d V — Iq, will be discussed in Sect. 9 below. 

The strong multiple absorptions seen in the Balmer 
lines and the Orion absorption system, so prominent at 
the time of second maximum, quickly vanished after the 
maximum. The only absorption lines still present in the 
December 10 spectrum are Nil 4601, NIII 4097-4103, OIII 
5592 A, and a weak, residual —1750 km/s absorption com- 
ponent in Ha. They too had vanished by December 12. 
The emission line spectrum during decline from second 
maximum changed much less than during decline from pri- 
mary maximum. The ionization/excitation increased only 
modestly, while there was a strong increase in equivalent 
width of the emission lines. The major change was seen in 
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Fig. 11. Spectroscopic evolution of Nova Cyg 2006 during second maximum in autumn 2006 (see Fig. 5 for plotting 
details). 



the optical depth of [OI] that soon started to decline from 
T6300~7 at maximum, to t 63 oo«0.15 on December 12. 

8.2. Decline after December 12 

Around December 12 the rapid decline of NCyg06 from 
the second maximum suddenly slowed down, settling to 
a <A£?>=0.0088 mag day" 1 straight linear slope over 
the range +300<A£<+460. This transition in the photo- 
metric evolution corresponds to the transition from stel- 
lar to nebular spectral types. In fact, the spectrum for 
December 10 in Fig. 11 is still dominated by the con- 
tinuum emission with a modest contribution from emis- 
sion lines. Three days later (cf. December 13 spectrum in 
Fig. 11), the spectrum had already turned into the nebular 
type, with a huge increase in the equivalent width of emis- 
sion lines and weakening of the underlying continuum. In 
particular, [OIII] 4363 was emerging on the red wing of 
the H7 profile, [OIII] 4959, 5007 were already key con- 
tributors to the strong emission blend with Nil centered 
at 4990 A, [Nil] 5755 was already stronger than H/3, and 
Hell 4686 was becoming detectable superimposed over the 
huge NIII emission at 4640 A. During the first weeks of 
the nebular state, nebular lines appeared double-peaked 
with a peak-to-peak velocity separation of 595 km/s for 
[OI] 6364, 815 km/s for [Nil] 5750, and 950 km/s for 
[OIII] 5007 A. 

With the transition from stellar to nebular type, 
the emission from the rapidly growing fraction of ion- 
ized ejecta compensated for the rapid photometric de- 



cline associated with the retracting pseudo-photosphere, 
which was increasing in temperature and therefore 
was shifting the peak of its emissivity to shorter and 
shorter UV wavelengths. The reprocessing into U BVRcIc 
bands and longer wavelengths of ultraviolet flux ab- 
sorbed by the ionized ejecta resulted in the leisurely 
<A_B>=0.0088 mag day -1 decline rate and in the pro- 
tracted visibility of the nova for such a long period. 
Without reprocessing ejecta, the central source would have 
become unobservable in the optical by January 2007. 



8.3. The fast and tenuous wind of mid 2007 



As illustrated by the spectra in Fig. 14, in July-September 
2007 the emission lines grew in width. The FWHM of 
[OI] lines increased from 1325 km/sec on March 11 spec- 
tra, to 1870 km/sec on July 20, and 1890 km/sec on 
September 19, and [Nil] lines from 1570 km/sec to 2070 
and 2135 km/sec, respectively. This trend is illustrated 
in Fig. 1, which also shows how the photometric decline 
almost stopped during the same period, with an extra 
~0.6 mag flux in the B and V bands with respect to ex- 
trapolated earlier lightcurves. We interpret these events 
as due to the emergence of a fast, tenuous wind from the 
central star. Its effects are equally evident in the evolution 
of the blended Ha+[NII] profile in Fig. 12. 
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Fig. 10. Middle panel: Photometric evolution of Nova Cyg 
2006 around second maximum. B-band lightcurve is at the 
top, Rc center, ic bottom. For telescope identification let- 
ters see sect. 2.1. Bottom panel: color evolution as recorded 
with telescope a. Upper panels: evolution of integrated flux 
and equivalent width of selected emission lines, expressed 
in magnitudes with respect to highest value. The 4660 
blend is initially dominated by NIII, than Hell become 
a relevant contributor. The 4990 blend is initially due to 
Nil, then [OIII] becomes the dominant contributor. 



9. Dust formation 

Infrared photometry of NCyg06 is available from the 
literature for only three dates: April 30, June 14, and 
December 12, 2006 (Russell et al. 2006; Mazuk et al. 2006; 
Lynch et al. 2006). These IR magnitudes were combined 
with our simultaneous UBVRcIc from Table 1 to con- 
struct the spectral energy distributions shown in Fig. 13, 
where the interpolating blackbodies have been reddened 
to the E B ^ V =0.56 affecting NCyg06. On April 30 and 
June 14, the UBVRcIc JHK spectral energy distribu- 
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Fig. 12. Evolution of Ha profile during the plateau phase, 
the second maximum and the advanced decline. At and 
AV are days and magnitudes from maximum. 



tion can be ascribed entirely to emission from the pseudo- 
photosphere with no contribution from warm dust, while 
on December 12 (At=+251 d ) thermal emission from dust 
is quite apparent. 

The dust detected on December 12 was optically thin. 
The energy re-emitted by dust in the infrared is energy 
absorbed from energy radiated by the nova, and therefore 
the dust optical depth is T^ImjL. L is difficult to es- 
timate directly at the time of dust condensation because 
most of it is radiated in the ultraviolet as indicated by the 
nebular and high ionization condition of the ejecta and 
illustrated by the spectral energy distribution of Fig. 13. 
However, theoretical models argue for a prolonged plateau 
period of constant luminosity experienced by novae fol- 
lowing optical maximum (Prialnik 1990, and references 
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therein), which is confirmed by multi- wavelength obser- 
vations of well studied novae (Gallagher 1977; Gallagher 
and Starrfield 1978; Gehrz et al. 1998). It can therefore be 
written that 

Tdust = t (6) 

-^plateau 

In Sect. 11 below, we show from photo- ionization model- 
ing that on January 20, 2007 (Ai=+290) the bolometric 
luminosity of NCyg06 was Mboi = —6.45 corresponding to 
Lbo\= 1.3xl0 38 erg sec -1 . We assume this was valid for 
December 12, 2006 as well. The luminosity radiated by 
the 1550 K blackbody fitting the December 12 dust emis- 
sion over the L, M, and TV-band photometric points of 
Fig.13 is L max = 2.6xl0 37 erg sec -1 . The optical depth of 
the dust for December 12 is therefore Tdust^O.2, indicating 
optically thin conditions. 

The infrared data available for NCyg06 do not enable 
a derivation of physical properties of the dust. However, 
for sake of discussion, we assume that the dust grains 
condensed in NCyg06 followed the mean properties ob- 
served in other novae (Gehrz 1988; Mason et al. 1996; 
Evans et al. 1997; Gehrz et al. 1998), i.e. they were small 
carbon grains (radius a<l /im, density p~2.3 gr cm -3 ), 
for which the Planck mean emission cross section goes as 
Q e = 0.01aTj MSt . Under these assumptions, the mass of 
the dust in NCyg06 on December 12 is 

M dust = 1.17 x 10Vr d -6 st (^j = 1.3 x 10 -9 M (7) 

which is quite small. Comparing with the amount of gas in 
the ejecta estimated below in Eq.(9), the gas-to-dust ratio 

in the ejecta on that date is ( ) ~ ^ x ^ sminar ly 
high value of the gas-to-dust ratio has been previously 
observed in Nova Vul 1984 N.2, which, like NCyg06, dis- 
played a large oxygen overabundance ([0/H]=+1.6, Gehrz 
et al. 1998). It should also be noted that dust did not ap- 
pear to form in some other novae, adequately observed in 
the infrared. 

As a first guess, it may seem that the hard radiation 
field from the central star, responsible for the ionization 
of the ejecta, should suppress dust grain formation (e.g. 
Gallagher 1997). Actually, observations of well-studied no- 
vae indicate that dust condenses when the ejecta become 
optically thin in the near ultraviolet and the nova spec- 
trum shifts from stellar to nebular type. Shore & Gehrz 
(2004) proposed that ionization may be a promoting agent 
of dust condensation via ionization-mediated kinematic 
agglomeration of atoms onto molecules and small grains 
through induced dipole interactions. It is interesting to 
note that dust was detected in NCyg06 on December 12, 
right when the transition was taking place. When compar- 
ing the December 10 and December 13 spectra in Fig. 11, 
the huge increase in both the amount and excitation de- 
gree of nebular emission is outstanding. Description of the 
infrared observations by Rayner et al. (2006) indicates 
that dust emission was absent on November 30 (when 
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Fig. 13. The optical and infrared spectral energy distribu- 
tion of Nova Cyg 2006 showing the presence of optically 
thin dust emission on December 2006. Optical data are 
from Table 1, and infrared data from Russell et al. (2006), 
Mazuk et al. (2006), and Lynch et al. (2006). The curves 
are fits with black bodies at the indicated temperatures 
and reddened by 75s_y=0.56 following a 7tV=3.1 stan- 
dard law. 



the nova was peaking at second maximum), and was still 
present on December 20 when the transition to a nebular 
spectrum was completed ([OIII] 4949 and [Nil] 5750 A 
appreciably stronger than H/3). They estimated a black- 
body temperature of 1410 K beyond 2 /im. Once corrected 
for the i?B_y=0.56 reddening affecting the nova, their es- 
timate is in fair agreement with our derived 1550 K for 
December 12 data. 

Such a temperature is unusually high for dust grains 
in novae (cf. Gehrz 2002). However, it has to be noted 
that on December 12 the dust had just begun to form in 
NCyg06 (a matter of a few days or even hours), it had 
a quite thin optical depth and it was quite low in mass. 
These are conditions away from the average ones applying 
to the novae for which dust properties have been studied in 
detail. We do not discuss the dust temperature in NCyg06 
further, since we lack the necessary infrared data. In par- 
ticular, contribution by emission silicate features at 10 /im 
and hydrocarbon features at 3.3, 3.4, and 11.3 /im have 
been observed in some novae superimposed onto the car- 
bon dust continuum during the optically thin phase (for 
ex. in V705 Cas by Gehrz et al. 1995, or V842 Cen by 
Gehrz et al. 1990). The presence of these (strong) emis- 
sion features may alter the slope of the spectral energy dis- 
tribution and consequently spoil the significance of tem- 
perature estimates from black-body fitting to broad-band 
photometric measurements. 

An intriguing feature of NCyg06 is the Ic band ex- 
cess it displayed during the second half of December (see 
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Fig. 10), soon after the dust detection on December 12. 
Following the second maximum, all bands declined steeply 
and in a similar fashion, with a marked slope change 
around December 12 when the transition from stellar to 
nebular spectrum took place and the ionization spread 
widely through the ejecta. From that date, all bands de- 
clined linearly with time, except ic, which displayed a pro- 
nounced extra-emission of A/c~l mag amplitude peaking 
around December 21 (see Fig. 10). During this time inter- 
val, we did not obtain spectra covering the Ic wavelength 
range. The closest spectra were taken on December 12, 
when the transition to nebular spectrum was already 
rapidly proceeding, and it shows only a modest 01 8446 A 
emission line and a weak Paschcn continuum in emission, 
neither of which are able to account for a sizable frac- 
tion of the Ic flux. We are therefore inclined to exclude 
a flaring of emission lines selectively concentrated over its 
wavelength range as the cause of extra emission in the Ic 
band during the second half of December. This is also sup- 
ported by the evolution of the V — Rc color illustrated in 
Fig. 10, which was dominated by the huge increase in the 
Ha equivalent width (cf. Fig. 1) induced by the nova tran- 
sition to the nebular stage (on the December 21 spectra, 
at the time of Ic band bump maximum, Ha was account- 
ing for 0.63 mag of the whole flux in the Rc band). The 
V — Rc color evolution is smooth and linear throughout 
the Ic bump development. 

The only feasible source of extra-emission that could 
have been responsible for the Ic December bump seems to 
be the dust. For dust emission to have leaked into the Ic 
band, it must have condensed into larger quantities and/or 
peaked to hotter temperatures than on December 12. A 
choice between the two alternatives could be possible only 
with the support of (missing) infrared data covering the 
second half of December. 

10. Comparison with V1493 Nova Aql 1999a 

The lightcurve of NCyg06 has been highly peculiar, but 
not completely unprecedented. A close match is repre- 
sented by Nova Aql 1999 N.l (=V1493 Aql), and Fig. 15 
compares the V-band lightcurves of the two novae. The 
similarity is striking, with V1493 Aql requiring 5.0 x less 
time to go from primary to secondary maximum, even if 
its t 2 ~7 and t 3 fa23 days (Venturini et al. 2004) arc not 
that different from t 2 =10A and i 3 =24 of NCyg06. Both 
novae seem to have orbital periods above the period gaps 
of cataclysmic variables, 5.0 hours for NCyg06 (Goranskij 
et al. 2006) and 3.7 hours for V1493 Aql (Dobrotka et 
al. 2006). Both novae belong to the Fell class and shared 
similar expansion velocities at first (FWHM=1700 km/s 
for V1493 Aql, Arkhipova et al. 2002) and second max- 
imum (FWHM=2000 km/s for V1493 Aql, Venturini et 
al. 2004). The spectral evolution of V1493 Aql has been 
poorly monitored. However, what is recorded suggests a 
path parallel to that of NCyg06: (a) low ionization con- 
ditions at both primary (only Balmer and Fell lines in 
emission in optical spectra, Tomov et al. 1999) and sec- 



ondary maxima (only Paschen, Bracket, OI and NI lines 
in emission with marginal Hel in IR spectra, Venturini 
et al. 2004), (6) slight increase in ionization between the 
two maxima with appearance of Oil, [Oil] , Nil and [Nil] 
emission lines and disappearance of Fell (Arkhipova et al. 
2002), (c) spectra obtained during the advanced decline 
characterized by nebular conditions, with [OIII], [Nil] and 
NIII among the strongest emission lines (Arkhipova et al. 
2002). 

The only possible major difference could be in the color 
evolution. Bonifacio et al. (2000) assembled heterogeneous 
B — V measurements published in several IAUC with some 
of their own, and concluded that V1493 Aql was hotter 
at maximum light (B — V<0.4) than at secondary peak 
(B — V^l.2) or plateau phase (B — V^0.9), with a bluer 
color during advanced decline (B — V~0.7). Unfortunately, 
the intrinsic variability of the comparison stars chosen 
by Bonifacio et al. to reduce their own measurements 
and the wild heterogeneity of the measurements collected 
from the literature does not call for any high confidence 
in the results. The reported color evolution could have 
been strongly affected by contamination to the V-band 
by the extremely strong Ha emission displayed by V1493 
Aql. Another word of caution about the Bonifacio et al. 
(2000) results is the excessively blue color they report for 
t 2 , B - V=0.32. After correcting it for the E B - V =0.6 de- 
rived by Venturini et al. (2004), the intrinsic color would 
be (B — V)o =— 0.3, equal to that of the hottest known 
stars {{B - V) o =-0.32 for 05 stars, Fitzgerald 1970), 
which contrasts with the cool temperature of the pseudo- 
photosphere required by the very low ionization condi- 
tion of the emission line spectrum. The discrepancy would 
be even greater if the i?s_y=1.5 reddening estimated by 
Arkhipova et al. (2002) would be adopted instead. The 
(B — V) =— 0.3 color at t 2 also sharply contrasts with the 
average (B — V) =— 0.02 derived by van den Bergh and 
Younger (1987) and also perfectly matched by NCyg06 (cf. 
Sect. 4 above). In conclusion, we are highly suspicious of 
the odd color evolution reported by Bonifacio et al. (2000) 
for V1493 Aql, and consider it worth further investigation. 

11. Photo-ionization analysis 

A photo-ionization analysis of NCyg06 spectrum for 
January 20, 2007 (Ai=+290 d ) has been performed with 
the CLOUDY code, version c06 (Ferland et al. 1998). This 
is the spectrum (plotted in Fig. 14) during the nebular 
phase with the highest S /N and lower line blending, before 
the flaring of the line profile widths seen at later phases 
(cf. Fig. 12). A detailed modeling of all collected spectra 
is far beyond the scope of this paper, and will be consid- 
ered elsewhere. Nevertheless, the results we obtained for 
the January 20, 2007 spectrum are confirmed by prelim- 
inary CLOUDY analysis of the other spectra, especially 
concerning chemical abundances. 

We assumed spherically expanding ejecta with a gas 
density declining as p(r)ocr~ 2 . We explored other density 
profiles (p(r)(xr n , with n=0, — 1, and —3), but obtained 
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Fig. 14. Spectroscopic evolution of Nova Cyg 2006 during the advanced decline from December 2006 to September 
2007 (see Fig. 5 for plotting details). 



a poorer match with observations. We did not fix inner 
and outer radii for the ionized shell (r in and r out ), and 
treated them as free parameters along with the ratio of 
the filled to the vacuum volumes in the ejecta (£). Only 
the abundances of chemical elements with observed lines 
were allowed to change, and all others were kept fixed to 
their solar value. Table 7 lists the observed and computed 
fluxes for the emission lines used in the modeling, Table 8 
summarizes the modeling results that are discussed be- 
low, and Table 9 provides the corresponding chemical mass 
fractions. 

11.1. Internal and external radii 

The shell of ionized gas at Ai=+290 d is a thin one, extend- 
ing from ri n —lA9 to r ou t=168 AU. The inner radius is a 
density boundary (no interior matter) . The outer one is an 
ionization boundary, defined by the distance at which the 
electronic temperature drops to negligible values (all ion- 
izing photons are absorbed within r out , and only neutral 
material could exist outward). These radii correspond to 
expansion velocities of 890 and 1010 km sec -1 . It is worth 
noticing that this is within the range of velocities displayed 
by the p\ principal systems of the multi-component ab- 
sorption spectrum, that it was the first to appear, the 
longest lived, and the one carrying the largest mass out- 



Table 7. De-reddened emission line ratios relative to H/3 
as observed on the Jan 20, 2007 spectrum of Nova Cyg 
2006 and as predicted by the CLOUDY model detailed in 
Table 8. 



line 


obs 


comp 


Ha + [Nil] 6548-84 


4.50 


3.40 


H/3 


1.00 


1.00 


H 7 + [OIII] 4363 


1.96 


3.27 


H<5 


0.36 


0.29 


He 




0.19 


H8 


0.20 


0.13 


5876 Hel 


0.14 


0.14 


7065 Hel 


0.06 


0.08 


4686 Hell 


0.20 


0.14 


4679 Nil 


0.06 


0.07 


5755 [Nil] 


1.95 


1.70 


4640 NIII 


0.38 




6300 [OI] 






7325 [Oil] 


0.70 


0.57 


4959 [OIII] 


1.32 


1.16 


5007 [OIII] 


4.30 


3.40 


3869 [Nelll] 


0.10 


0.12 


4720 [NelV] 


0.005 


0.0045 


7135 [Arlll] 


0.016 


0.016 


7236 [ArlV] 


0.0013 


0.0016 


7006 [ArV] 


0.0024 


0.0023 


6087 [FeVII] 


0.007 


0.007 
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Fig. 15. Comparison of Nova Cyg 2006 F-band lightcurve 
with the m vis lightcurve of Nova Aql 1999 N.l (AAVSO 
visual estimates). The ordinate scale is the same in both 
panels. The time scale of Nova Aql 1999 N.l abscissas is 
5.0 x faster that that of Nova Cyg 2006. 

Table 8. Parameters of the CLOUDY model best-fitting 
the de-reddened emission line ratio of Table 7. From top to 
bottom: temperature and radius of the central black-body 
source, inner and outer radii of the emitting shell, hydro- 
gen density at inner and outer radii, electronic density at 
the inner and outer radii, ratio of the filled to vacuum 
volumes, clement abundances relative to solar. 



log 21| (K) 


5.236 


log Ebb (cm) 


10.14 


logr in (cm) 


15.35 


logr out (cm) 


15.40 


\ogpg (cm" 3 ) 


7.781 


lo SP™t (cm" 3 ) 


7.874 


logiVj^ (cm -3 ) 


8.00 


log^t (cm" 3 ) 


5.00 




0.12 


H/H 


0.80 


He/He 


0.87 


N/N 


316 


o/o 


26 


Fe/Fe 


1.0 


Ar/Ar 


1.0 


Ne/Ne 


1.0 



Table 9. Mass-fraction abundances of measured elements 
in Nova Cyg 2006 and, for reference, in the Sun. 



mass-fractions 



X 

Y 

Z 

N 

O 

Ne 

Ar 

Fe 



NCyg 

0.427 

0.185 

0.388 

0.219 

0.163 

0.0013 

0.00008 

0.00096 



Sun 

0.704 

0.279 

0.017 

0.001 

0.008 

0.0017 

0.00011 

0.00126 



flow, as Fig. 6 and Table 4 clearly illustrate. Therefore, 
the ejected material causing the p\ principal absorption 
systems during early post-maximum decline appears to 
be the same, creating most of the emission line fluxes at 
At=+290 d during the nebular phase. 

11.2. Mass in the shell 

The hydrogen mass fraction in the ionized shell is 
A=0.427, and the ratio of the filled to vacuum volumes is 
£=0.12. Therefore, the total gas mass within the ionized 
shell at At= 



=+290 d is 



M. 



shell 



X h in 



4irr 2 p(r)dr = 3 x 10~ 4 M c 



© 



(8) 



In Sect. 5.3.3 we estimated from the [OI] lines the amount 
of neutral gas in the ejecta during the early days of the de- 
cline from maximum when the ejecta were predominantly 
neutral, and found M^ cd = 5 xl0~ 4 M Q (Eq. 5). The 
agreement of the two estimates is quite reassuring, in par- 
ticular if it is taken into account that the ionized shell 
at At=+290 d is radiation-bounded and neutral material 
could exist outward. Therefore M s hcii = 3 xl0~ 4 M© is a 
lower limit to the total amount (ionized + neutral) of gas 
existing at At—+290 d . For sake of discussion in the rest 
of the paper, we assume for the total amount of ejected 
material an averaged value between the two methods: 



<cc a tcd = 4 x 10- 



11.3. Central star 



Mr, 



(9) 



The central ionizing source is found to have a radius 
#=0.198 R , a temperature T eff =172 000 K, and therefore 
a luminosity 3.08xl0 4 L Q , corresponding to Mb i=- 6.45 
(assuming Mboi=4.74 for the Sun, Livingston 2000). 

11.4. Mass of the white dwarf 



At the time of the January 20, 2007 spectrum, NCyg06 
was undergoing the luminosity plateau phase. It char- 
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acterizes the observed and theoretical evolution of no- 
vae after the initial super-Eddington phase and before 
the switching-off of the nuclear reactions and consequent 
cooling track for the white dwarf (cf. Gehrz et al. 1998). 
During the luminosity plateau phase, the white dwarf 
burning residual hydrogen at the surface follows the 
Paczyhski (1971) core-mass-luminosity relation L cm \ = 
6 x 10 4 (M W d/M q - 0.522) L©. Inserting the above es- 
timated luminosity, the mass of the white dwarf turns out 
to be 

M WD = 1-03 M Q (10) 

The amount of material to be accreted on this white dwarf 
for the pressure to reach 10 20 dyne and trigger a TNR 
event is 4xl0~ 4 M Q (cf. Starrfield 1989). In comparison 
with Eq.(9), the outburst therefore seems to have ejected 
most of material accreted by the white dwarf before the 
eruption set in. 

11.5. Chemical abundances 

The chemical abundances in the ionized shell are given in 
Table 8 (relative to solar) and Table 9 (as mass fractions) . 
They clearly reflect the non-equilibrium CNO-cycle burn- 
ing of hydrogen, with massive overabundance of oxygen 
and nitrogen, among the largest ever measured in novae. 
Carbon should be enhanced too, but there were no de- 
tectable carbon emission lines in the January 20, 2007 
spectrum of NCyg06 that could be used to measure its 
abundance. 

Iron is not enhanced. Because it is not produced by 
the nuclear reactions, it is expected to reflect the metal- 
licity of the donor star, which should be around solar given 
the likely association of NCyg06 with the galactic disk at 
72=11.4 kpc galacto-centric distance. 

The large overabundances of Ne observed in the spec- 
tra of some of the fastest novae are ascribed to mixing 
into the accreted envelope of material from the underly- 
ing white dwarf (e.g. Gil-Pons et al. 2003). In massive 
progenitors, non-degenerate carbon ignition leads to the 
formation of a degenerate core mainly made of oxygen and 
neon. The minimum mass on the zero age main sequence 
leading to extensive carbon-burning is M~9.3 M and the 
resulting white dwarf will have a mass of Mwd>1-1 M q . 
Neon is not enhanced in NCyg06, in agreement with the 
above estimate of the white dwarf mass, which is below 
the limit for the ONe type. 

12. Energy budget 

The total energy budget of the outburst is the sum of the 
energy radiated from the central star and the binding and 
kinetic energy of the ejecta. 

The energy radiated by the central star was computed 
by integrating the bolometric luminosity over time. For a 
rough order of magnitude estimate, we assumed that the 
Mboi=— 8.0 value at maximum declined linearly with time 
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to Mboi=— 6.45 on At=+290 d , remaining there until the 
end of our observing campaign on At=+530 d . The result 
is 

,•530 

Erad = / M bol {t)dt = l x 10 46 ergs (11) 
Jo 

The actual time evolution of the bolometric luminos- 
ity is hard to define. Before NCyg06 settled to the 
constant luminosity phase, represented by Mb i=- 6.45 
on At=+290 d , it went through a short-lived, super- 
Eddington phase characterized by Mboi=— 8.0 at peak V- 
band brightness. This is a common feature of observed no- 
vae and theoretical models (e.g. Starrfield 1989). Contrary 
to most novae, NCyg06 went through a second maximum, 
and this could have been reflected in an increase in Mbol 
above its plateau level. Thus a linear trend in luminosity 
evolution is adequate to derive an order of magnitude for 
the radiated energy, which is moreover a marginal term in 
the overall energy budget. 

The kinetic energy of the ejected shell is 

E kin = ^ jf 2nr 2 p{r) dr = 2.2 x 10 47 ergs(12) 

where p{r) — p(r- m )(r /r- m )~ 2 and the velocity of the ejecta 
is taken to be v(r) = r/t which is consistent with an ejec- 
tion over a short period of time around t=0. 

The outburst also had to provide the energy required 
to gravitationally unbound the ejected material. The bind- 
ing energy depends from the mass of the white dwarf (its 
potential well). The binding energy of the above estimated 
M e t ? t c a t 1 cd =4 x 10~ 4 M from a 1.03 M white dwarf is 

E bm = 2 x 10 47 ergs (13) 

which is quite similar to the kinetic energy. As is typical of 
novae, most of the outburst energy goes into the mechan- 
ical work of mass ejection, and only a few percent makes 
up the radiated fraction. 

The total energy released by the outburst is therefore 

E tot = E bm + E rad + E Hn = 4.3 x 10 47 ergs (14) 

which corresponds to the hydrogen burning of 
4.6xl0~ 5 M Q of accreted material of solar composi- 
tion, which is about 10% of the amount of mass accreted 
by the white dwarf prior to the outburst. 
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